Infinitely Robust Order and Local Order-Parameter Tulips 
in Apollonian Networks with Quenched Disorder 
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For a variety of quenched random spin systems on an Apollonian network, including ferromagnetic 
and antiferromagnetic bond percolation and the Ising spin glass, we find the persistence of ordered 
phases up to infinite temperature over the entire range of disorder. We develop a renormalizat ion- 
group technique that yields highly detailed information, including the exact distributions of local 
magnetizations and local spin-glass order parameters, which turn out to exhibit, as function of 
temperature, complex and distinctive tulip patterns. 



Although their structure dates back to ancient Greek 
mathematics, Apollonian networks |2] have seen a 
recent surge of interest as a simple and elegant model 
that incorporates some of the key features identified in 
real-world networks: a scale-free degree distribution, the 
small- world effect, and a high clustering coefficient. As 
such, they have become a versatile tool for understanding 
the effects of complex topologies in interacting systems: 
applications include percolation and epidemic spread- 
ing [l|, 0], mametic systems Q, mechanisms of net- 
work growth J5| , avalanches in sandpile models [6] , neu- 
ral networks [3] , and even quantum behaviors like coher- 
ent exciton transport [8] and correlated electron mod- 
els jol, [13] . The latter were inspired by the development of 
synthetic, nanoscale, non-branched fractal polymers 
which raises the possibility that the unique properties 
of scale-free structures like Apollonian networks may be 
harnessed for technological applications. 

In this study we focus on an intriguing aspect of 
these networks: their ordering resilience in the pres- 
ence of imposed quenched disorder. For an Ising model 
with a variety of random-bond distributions - ferro- 
magnetic/antiferromagnetic percolation and spin-glass - 
we find that ordered phases persist up to infinite tem- 
perature, for every case except the pure antiferromag- 
netic system (where the geometrical frustration of the 
network leads to paramagnetism). The self-similar na- 
ture of the Apollonian network allows the use of ex- 
act renormalization-group (RG) techniques to calculate 
the phase diagram structure, even in the presence of 
quenched randomness. While there have been many nu- 
merical RG studies of spin glasses on fractal lattices, 
we have gone further, developing an iterative procedure 
based on the local recursion matrix that allows us to 
calculate exactly the complete distribution of the local 
magnetization and spin-glass order parameters, for the 
full range of bond probabilities and temperatures. The 
resulting local-order diagrams (Fig. [1]) show an intricate 
structure as temperature is lowered, never before ob- 
served in such detail for a disordered spin system. Apol- 



lonian networks can be embedded in a Euclidean plane 
without any edge crossings pj, l2| . This makes spin sys- 
tems on such networks potentially physically realizable, 
for example in nanostructures formed from dense poly- 
disperse packings of magnetic grains [1, 4]. Apollonian 
packings have also been used in the study of smectic liq- 
uid crystals [l2[. 

Network structure - The construction of the Apollonian 
network is depicted in Fig. [2l^a): at each step, a new site 
is added to the center of every triangle in the network, 
and connected to the surrounding vertices. In the limit 
of infinite size, the geometrical characteristics of the net- 
work can be summarized as follows 

0,0,0: P{k) being 
the probability that a site has degree /c, the cumulative 
degree distribution is Pcum(^) = YlV'=k-P(^) ^ 
large /c, with the scale- free exponent 7 = 1 + In 3/ In 2?^ 
2.585. Due to the compact network structure, the aver- 
age shortest-path length I between any two points scales 
as in the small- world effect, I ~ In A/, as shown in 
using the exact recursive method of [14] . As is typical in 
small- world networks, the average clustering coefficient is 
large, C ^ 0.828, measuring the ratio of the connections 
among the nearest-neighbors of a site and the maximum 
possible number of such connections /c(/c + l)/2, where k 
is the degree of the site. In all these respects the topolog- 
ical properties of the Apollonian network are comparable 
to those observed in empirical complex networks that are 
simultaneously scale- free and small- world 

We study an Ising Hamiltonian on the network, 
-PH = JijSiSj, where Si = ±1, P = l/ksT, 

the sum {ij) is over nearest neighbors, and the bond 
strengths Jij are distributed with a quenched random 
probability distribution P{Jij). We consider two types 
of distributions described by bond probability p : (i) 
the percolation case, where P{Jij) = p6{Jij — J) -\- {1 — 
p)S{Jij)^ for both ferromagnetic (F) J > and antiferro- 
magnetic (AF) J < interactions; (ii) the ± J spin glass 
(SG) case, where P{Jij) = pS{Jij + J) + (1 -p)S{Jij - J) 
with J > 0. 

Exact renormalization-group transformation - The self- 
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FIG. 1: (Color) Local order-parameter tulips: Probability distributions of local magnetization (left panels) and spin-glass 
(right panels) order parameters of the interior sites on an Apollonian network with Ising spin-glass interactions, as a function 
of temperature, for three different antiferromagnetic bond concentrations p. 

are functions of the original interactions Ji , J2 , J3 of the 

.^.^ 7^ _ l7 [ cOsh(2Ji+2J2)+COSh(2J3) l . y 

Star, J12 - -^in | ^^^^ (2Ji-2J2)+cosh (2J3) / ^^^^^^ 
permutations. In the context of the whole network, the 
mapping works as shown in Fig. EJc). The Jij interac- 
tions (inner to each triangle) are added to the original 
interactions Jij of the (n — l)th generation network (in- 
between triangles), to give the renormalized interactions 
J^j. Thus, in the bulk each original interaction gets Jij 
contributions from its two adjoining plaquettes. 

In order to implement this RG transformation for the 
system in the thermodynamic limit, we focus on the prob- 
ability distribution of triplets Q{{Jij, Jjk, Jik}) gener- 
ated by the star-triangle transformation. As we iterate 
the RG mapping, this distribution Q changes, and we 
can extract thermodynamic information from the result- 
ing flows. To keep track of Q at each step, we adapt a nu- 
merical procedure developed by Nobre [l7| for RG trans- 
formations of spin glasses on hierarchical lattices. This 
method has been shown to give numerically accurate re- 
sults for phase diagrams [18], agreeing with more compli- 
cated binning techniques used to directly evaluate the RG 
flows of interaction distributions [13]. We represent the 
distribution Q by a pool of large size M, where each ele- 
ment in the pool consists of a triplet of real numbers. To 
generate the initial pool we repeat the following M 
times: (i) choose three random numbers Ji, J2,-^3 with 
the probability P( J); (ii) perform the star-triangle trans- 
formation of Fig. [2jb), yielding a triplet { J12, ^23, ^13} 
which is placed in the pool. Each subsequent RG trans- 
formation creates a new pool Q^*^ from the previous pool 
in the following manner, again repeating the same 
procedure M times to preserve the size of the pool: (i) 




FIG. 2: (Color online) (a) Construction of an Apollonian net- 
work, (b) Star-triangle transformation, (c) Two successive 
RG transformations of an Apollonian network. 



similar structure of the Apollonian network allows us to 
formulate an exact RG transformation. For —pl-t{{Jij}) 
the Hamiltonian for a particular configuration of interac- 
tions on the nth-generation network (i.e., the lattice after 
n construction steps), the RG mapping yields a Hamil- 
tonian —P'H' {{J'ij}) with a renormalized set of interac- 
tions {J'ij} on the (n — l)th-generation network, preserv- 
ing the partition function. The mapping is carried out 
through a star-triangle transformation [16|, tracing over 
the spins at sites added at the nth step. This is shown 
for one plaquette in Fig. [2](b), with the decimated spin 
labeled a. The trace over <j produces interactions J12, 
-^13, -^23, between the edge sites of the triangle, which 
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FIG. 3: (Color online) Phase diagram of the Ising spin glass 
on an Apollonian network, in temperature 1/J versus antifer- 
romagnetic bond concentration p. The boundary between the 
ferromagnetic and spin-glass phases is first order. The para- 
magnetic phase appears with a first-order phase transition at 
p^l. 



randomly choose three triplets from Q^*"^^; (ii) randomly 
arrange these three triplets like the three triangles in the 
second step of Fig. [2](c), together with the three mid- 
dle bonds, chosen randomly with probability P{J); (iii) 
decimate the center spin to yield a renormalized triplet, 
namely the inner bonds in the third step of Fig. [21(c)), 
which is placed in pool Q*^*^. As M ^ 00, the pools 
mimic the exact renormalized distributions of triplets 
in the thermodynamic limit. For the present work, we 
found that M = 10^ was sufficiently large to make finite- 
ensemble effects negligible. From the behaviors of the 
Q^*^ in the limit of large we can identify the phase 
structure of the system. Specifically, looking at the av- 
erage J*^*^ and standard deviation ctj^ of the 3M bond 
strengths in pool we can distinguish three limiting 
behaviors as i ^ 00: (i) a ferromagnetic (F) sink, where 
J*^*^ 00, cTj^ 00, G^j' / J^^^ 0; (ii) a spin-glass 



(SG) sink, where J^^ 



00, (J 



00 



0; 

(iii) a paramagnetic (P) sink, where J^*^ 0, ctj^ 0. 
Furthermore, the RG evolution of fraction of frustrated 
triangles is different in the ferromagnetic and spin-glass 
phases, respectively going to and 0.5 in the ferromag- 
netic and spin-glass phases. 

Calculation of local magnetizations and local SG order 
parameters - Moreover, the numerical procedure de- 
scribed above is not limited to just the distribution of 
renormalized interactions Q. It can be extended to deter- 
mine additional thermodynamic details, in particular the 
distribution of local magnetizations and local SG order 
parameters. Let us consider the magnetization rria- at a 
site a in the original lattice. For simplicity, let a be one of 
the sites generated at the last construction step. We shall 
denote these as "interior sites", and they constitute 2/3 
of the total lattice in the limit of large n. Adding a local 
magnetic field i^cr, we can write rria- = dlnZ/dHcj\H^=o^ 
where Z is the partition function. The star-triangle 
transformation with Her produces additional interactions 
in the renormalized triangle on the right side of Fig.[2](b): 
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FIG. 4: (Color online) Magnetization m and spin-glass or- 
der parameter q for the Ising spin glass on an Apollonian 
network. In the left panels, the curves, consecutively left- 
most from the top, are for constant temperature 1/J = 
0.2,1.0,2.0,10.0. In the right panels, the curves, consecu- 
tively from the top, are for antiferromagnetic bond concen- 
tration p = 0.10, 0.38, 0.50, 0.90. 



three local fields iJi^i, iJ2S2, ^3^3, and a three-site term 
KS1S2SS, where and K are functions of the Jij and 
Her. With these interactions, the RG mapping is closed 
upon further iteration, so there will be a set of parame- 
ters K^^^ = {H^\h!^\h^\k^^} after the iih RG step 
associated with the triangle that originally contained spin 
<j. Using the chain rule, rUcr can be expressed [l9|, [20| in 
terms of local recursion matrices T over the i steps. 



rUc 



where /ii^^ = d\iiZ/dK^J;\ T^l = dxf / 

dKa'' /dH(j^ and Ka'' are the components of K*^*^. All 
these quantities are evaluated in the RG subspace with 
initial condition iJ^r = (i.e., where Ka^ = for all 
z), which makes them functions only of the Jij configura- 
tion at the previous step. In the thermodynamic limit, as 
a corner boundary condition, we calculate /l^^*^ over the 
up-magnetized configurations of the three original corner 
spins of the network. As the RG transform is iterated, 
each triplet in the pool Q^^^ has a corresponding vector 
T(^) . . .t(2)v(^), which can be contracted with fi^'^ cal- 
culated from Q*^*^ to obtain a pool of rria- using Eq. ([1]). 
For sufficiently large i and M, the resulting pool con- 
verges to the exact distribution of local magnetizations 
in the thermodynamic limit. The averages of m^r and 
Qa- = over this distribution respectively yield the mag- 
netization m and SG order parameter q for the interior 
sites. 

Results - We focus first on the ferromagnetic (F) and 
antiferromagnetic (AF) percolation cases. For F percola- 
tion the system is ferromagnetically ordered at all finite 
temperatures for any p > 0, directly related to the ex- 
istence of a giant connected component in the network 
at all nonzero bond probabilities [i]. In contrast, one 
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might expect macroscopic order in the system to be in- 
hibited in the AF case, since AF bonds are frustrated 
on the triangular plaquettes in the network. Indeed, in 
the case of p = 1, namely for a pure AF system, frustra- 
tion leads to a paramagnetic phase at all temperatures. 
However, as soon as even a tiny fraction of AF bonds 
is removed, namely for any < p < 1, we find an SG 
phase at all temperatures^ an interesting example of a 
glassy phase which is completely impervious to thermal 
excitations. The SG phase appears even when weaker 
forms of disorder are added to the pure AF system, such 
as simply attenuating a fraction of the bonds. Consider 
the range of models described by the bond distribution 
P{Jij) = pS{Jij — J) + (1 —p)S{Jij —cJ) where J < and 
< c < 1. Here c = corresponds to the AF percolation 
case described above, but it turns out that any c < 1 
gives the same phase diagram structure: a SG phase of 
infinite extent for all < p < 1 and paramagnetism at 
p = l. 

We now turn to the spin-glass system: The system 
composed of antiferromagnetic bonds, under infinitesi- 
mal doping by ferromagnetic bonds, produces a spin-glass 
phase via a jump, namely a first-order phase transition 
at p = 1. For a sufficient quantity of ferromagnetic dop- 
ing, a first-order transition occurs to the ferromagnetic 
phase, as can seen in the phase diagram in Fig. [31 On 
both sides of the transition line the ordered phases per- 
sist to infinite temperature 1/ J, and the boundary itself 
asymptotically approaches a vertical line at p = 0.5 as 
1/J oo. The first-order nature of the ferromagnetic- 
spinglass phase transition, to our knowledge not seen in 
other systems, is evident from the magnetization and SG 
order parameter plotted in Fig.Hl which indeed show dis- 
continuities crossing the boundary. (At the highest tem- 
perature depicted, 1/J = 10.0, the discontinuities exist 
but are too small to be seen on the scale of the figure.) 
Fig. m also reveals a curious aspect of spin-glass order 
on the Apollonian network: unlike a conventional spin- 
glass phase, the magnetization m is generally nonzero. 
With the above-mentioned boundary condition on the 
three corner spins, m ^ for p ^ 0.5. The negative m 
at large p is understood from the infiuence of the corner 
spins, which except for the central spin have the highest 
degree in the network. In an environment with mostly 
antiferromagnetic bonds, an up orientation for the corner 
spins will yield a negative magnetization of the interior 
sites. This ability of the most connected spins to de- 
termine the sign of the magnetization may be a general 
feature of scale- free networks, and has been seen in the 
Barabasi- Albert model [2l[. 

The evolution of the system under disorder is obtained 
in microscopic detail by the calculation of the local mag- 
netizations and local SG order parameters, as described 
above. The resulting full distributions of the local magne- 



tizations and SG order parameters are given in Fig.[TJ To 
produce these graphs, the m^r and Qo- pools were coarse- 
grained using a binning procedure, and the normalized 
heights of the resulting histograms color-coded. For clar- 
ity, histograms smaller than 10~^ are not shown. The 
distributions exhibit a distinctive tuliplike shape, devel- 
oping a rich structure as the system is cooled, spreading 
from narrowly localized peaks at high temperatures into 
complex bands of smaller peaks over the whole range at 
intermediate temperatures. These bands in turn con- 
verge toward the expected sharply defined values at low 
temperatures, with local magnetization peaked around 1, 
0, and -1, and the SG order parameter around and 1. 
The asymmetry leading to negative magnetization m for 
p > 0.5 is evident in comparing the p = 0.5 and p = 0.9 
local magnetization plots. The former is entirely sym- 
metric between negative and positive peaks, while in the 
latter the predominance of antiferromagnetic bonds leads 
to the bands of negative peaks becoming more prominent. 

In conclusion, we have shown that ferromagnetic 
phases and, moreover, spin-glass phases on Apollonian 
networks exhibit a remarkable robustness, with an infi- 
nite critical temperature for any amount of disorder. In 
fact, order persists to infinite temperature even when al- 
most all of the bonds in the system are removed and even 
when almost all of the bonds in the system are frustrated, 
namely in the p infinitesimally greater than zero and in 
the p infinitesimally less than one regimes of the perco- 
lation and spin-glass problems, respectively. This prop- 
erty should have consequences for actual applications on 
networks. For example, interacting objects arranged on 
Apollonian nanostructures would be able to maintain co- 
operative behavior over a broad range of temperatures 
and intrinsic disorder. Our local renormalization-group 
theory method yields, in the network with frozen disor- 
der, the exact local order parameters, up to now only 
calculated approximately by mean-field theory. The re- 
sulting local magnetizations and local spin-glass order 
parameters do not yield just a distribution of values, as 
would be most simply expected in systems with frozen 
disorder, but also unexpectedly distinctive tulip struc- 
tures with stalks, leaves, and veins, evolving under tem- 
perature. It is clearly unlikely that the new and intrigu- 
ing tulip structures, with stalks, leaves, and veins, in 
the microscopies are limited to Apollonian networks, but 
more likely will appear, perhaps in varying topologies, 
in diverse small- world systems. Our locally discriminat- 
ing renormalization-group technique, in yielding such de- 
tailed local results, should be of interest for the positional 
distribution of order in systems with inhomogeneities, be 
it due to quenched impurities or surfaces, etc. 
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